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Smooth  Sliding  Mode  Controller  Design  for  Robust  Missile  Autopilot 


Description  of  the  Effort 

a.  General 

The  University  of  Alabama  in  Huntsville  (UAH)  will  provide  the  personnel  and 
equipment  to  develop  Smooth  Sliding  Mode  Control  algorithms.  These  algorithms  may  be  used 
to  provide  a  generic  flight  control  to  advanced  guided  weapon  systems. 

b.  Requirements 

The  developed  Sliding  Mode  Controller  (SMC)  algorithms  will  be  robust  to  agile  reference 
commands  and  unknown  disturbances  from  environmental  conditions  (atmospheric  turbulence), 
as  well  as  missile  and  its  actuator  uncertainties  and  failures.  The  key  task  in  developing  a  smooth 
SMC  is  to  design  a  multi-loop  control  system,  where  the  sliding  mode  exists  in  every  loop.  A 
virtual  control  signal  to  be  produced  in  every  outer  loop  must  be  smooth,  for  it  has  to  be  tracked 
by  inner  cascades  of  a  multi-loop  system.  Smoothness  of  a  control  signal  and  existence  of  a 
multi-loop  SMC  must  be  achieved  without  degrading  high  precision  and  robustness  of  sliding 
modes.  The  developed  smooth  sliding  mode  control  algorithms  are  to  be  applied  to  generic 
missile  control  problems. 
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Abstract 


Presented  is  a  method  of  smooth  sliding  mode  control  design  to  provide  for  the  second- 
order  sliding  mode  on  the  selected  sliding  surface.  The  control  law  is  a  nonlinear  dynamic 
feedback  that  in  absence  of  unknown  disturbances  provides  for  finite-time  convergence  of  the 
second-order  reaching  phase  dynamics.  The  application  of  the  second-order  disturbance  observer 
in  a  combination  with  the  proposed  continuous  control  law  gives  the  second-order  sliding 
accuracy  in  presence  of  unknown  disturbances  and  the  discrete-time  control  update.  The 
piecewise  constant  control  feedback  is  smooth. 

The  proposed  control  scheme  is  applied  to  a  generic  three-loop  control  system,  and  the 
proposed  multiple  time  scale  sliding  mode  design  is  demonstrated  on  the  simplified  model  of  a 
ballistic  interceptor  missile. 
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1.  Introduction 

The  area  of  advanced  missile  guidance  and  control  is  a  promising  field  of  application  for 
latest  developments  in  robust  control  theory.  Many  generations  of  highly  successful  missiles  use 
different  variants  of  classical  algorithms  such  as  proportional  navigation  guidance  and  PID 
control  [1-3]  in  association  with  various  extended  Kalman  filtering  techniques.  This  is  the  major 
reason  for  industry  designers  to  continue  advancement  and  modernization  of  these  relatively 
simple  topologies  and  linear  autopilots  with  fixed  given  structure.  However,  “emerging  threats 
and  new  operational  constraints  have  created  the  need  for  a  variety  of  new  weapon  systems, 
even  existing  systems  can  expect  significant  upgrades. ”[4]  The  demand  for  engaging  the  highly 
agile  targets  [5]  with  hit-to-kill  accuracy  and  resistance  to  electronic  countermeasures,  combined 
with  that  for  increased  standoff  range  or  high  off-boresight  angle  and  effective  reconfiguration 
for  different  missions,  create  new  challenges  for  future  missile  design.  There  one  can  observe 
three  main  directions  of  rapid  theory  and  technology  proliferation  to  address  the  need: 

•  advanced  navigation  equipment  and  estimation  algorithms; 

•  advanced  guidance  strategies; 

•  advanced  missile  airframes. 

In  the  last  direction,  the  high  performance  vehicles  are  designed  to  be  open-loop  unstable,  to  be 
non-axisymmetric  for  extreme  agility;  new  control  actuating  mechanisms,  including  hot-  or  cold- 
gas  thrusters  and  thrust  vectoring,  hybrid  rocket  motors,  reduced  surface  area  high-speed  fins 
etc.,  are  in  development  to  increase  maneuverability  over  expanded  envelope  of  flight 
conditions. 

In  many  cases,  integration  of  new  guidance  strategy  with  new  airframe  (hybrid  rocket 
motor  with  LQG  optimal  guidance  law  [6],  Sliding  Mode  Control  based  guidance  with  pulse 
type  on-off  thrusters  [7])  directly  accounts  for  increased  performance  in  missile  airframe  due  to 
the  mentioned  above  developments  only,  leaving  the  autopilot  structure  to  be  the  most 
conservative  area  for  design  upgrades,  and  frequently  reducing  the  tracking  control  problem  to 
control  allocation  only  [6]. 

Missile  autopilots  are  usually  designed  to  track  acceleration  issued  by  a  guidance  law  on 
the  basis  of  data  sensed  and  estimated  by  a  navigation  system,  where  a  widely  accepted  approach 
is  to  estimate  relative  position,  velocity  and  acceleration  via  extended  Kalman  filtering 
technique.  The  typical  autopilot  uses  the  gain  scheduling  technique  or  rather  different  parameter 
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optimization  methods  (see  a  review  in  [4]),  given  a  particular  fixed  controller  structure,  such  as  a 
classical  PID  architecture.  The  simplest  control  structure  is  easy  to  implement,  to  time  up  and 
validate  for  different  performance  specifications.  The  industry  designers  prefer  it  to  recent  linear 
and  nonlinear  schemes,  because  they  never  really  know  how  these  novel  techniques  of  high 
complexity  or  high-order  dynamics  would  trade  off  performance  for  uncertainty  robustness,  or 
improve  them  both  significantly,  or  whether  they  will  work  at  all.  Moreover,  many  design 
approaches  pursue  the  idea  that  the  problem  to  design  a  good  controller  is  to  obtain  a  good 
mathematical  model  of  the  plant,  to  derive  sophisticated  equations,  which  can  include  as  much 
known  effects  as  possible  and  utilize  as  much  as  possible  the  data  obtained  via  hardware-in-the- 
loop  simulations. 

Having  said  that,  one  can  come  to  the  following  conclusion. 

1.1  Requirements  for  Advanced  Missile  Autopilot 

A  new  missile  autopilot  of  practical  interest  should  have 

•  A  simple  clear  structure,  which  can  be  readily  designed,  facilitated,  implemented,  tuned  up 
and  monitored; 

•  Increased  both  performance  and  robustness  to  visible  level; 

•  Required  the  least  possible  modeling  information  to  guarantee  these  performance  and 
robustness; 

•  Relaxed  the  sensor/processor/actuator  requirements  without  introducing  new  restrictive 
assumptions,  which  is  difficult  to  verify; 

•  Maximal  autonomy  with  respect  to  common  sensed  data  to  be  easily  integrated  with  different 
navigation/guidance  systems; 

•  Required  all  on-the-fly  ad  hoc  inevitable  modifications  being  conducted  by  a  professional 
with  common  relevant  control  theory  background. 

Having  almost  all  the  mentioned  above  characteristics,  the  Sliding  Mode  Control  (SMC) 
systems  [8-12]  are  valued  in  the  class  of  robust,  high  fidelity  control  algorithms  for  their  robust 
accommodation  of  uncertainties  and  cancellation  of  disturbances.  In  addition,  a  new  generation 
of  sliding  mode  algorithms  (high  order  sliding  modes  [13], [14])  featuring  usually  1st  or  2nd  order 
dynamics  in  traditionally  static  nonlinear  state  feedback  (conventional  sliding  mode),  produces 
continuous  control  input,  which  can  explicitly  account  for  not  only  actuator  saturation  but  also 
rate  saturation  limit.  The  other  additional  features  include  enhanced  robustness  to  measurement 
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noise  and  opportunity  to  realize  output  feedback  instead  of  state  feedback  to  increase  robustness 
to  parametric  and  dynamic  uncertainties. 

1.2  SMC  as  a  Perspective  Control  Technique  for  Advanced  Missile  Autopilot 

Further,  we  briefly  outline  the  SMC  design  specifications  to  address  some  of  the 
formulated  requirements.  A  SMC  system  has  a  simple  controller  structure,  clear  distinct  design 
steps  and  performance  assessment.  It  starts  with  a  highly  uncertain  plant  model,  where  the 
robustness  analysis  (Lyapunov  function  technique)  is  included  in  controller  design.  The  only 
known  information  required  is  the  relative  degree  of  each  commanded  output  and  the  ranges  of 
uncertainties  for  unknown  parameters,  nonlinear  terms,  external  disturbances  and  measurement 
noise,  as  well  as  actuator  limits  and  rate  limits.  Uncertain  terms  in  certain  subspace  can  be  non¬ 
smooth.  For  noisy  signals,  the  value  of  the  second  time  derivative  above  which  a  signal  can  be 
considered  as  noise  is  required. 

The  first  step  in  controller  design  is  to  judiciously  select  the  sliding  manifold  (sliding 
surface)  in  the  system  state  space.  Equations  of  sliding  surface  and  the  internal  dynamics  of  the 
plant  compose  the  closed-loop  (compensated)  dynamics  of  the  plant.  In  minimum  phase  case, 
equations  of  sliding  manifold  are  the  desired  tracking  error  dynamics,  which  is  uncoupled  from 
the  internal  dynamics  and  completely  plant  independent.  The  system  motion  on  this  manifold 
(the  sliding  mode)  will  be  absolutely  robust  to  plant  uncertainties  and  external  disturbances.  The 
order  of  asymptotic  or  finite-time-convergence  tracking-error  dynamics,  which  can  be  enforced 
in  sliding  mode  is  equal  to  r-1,  where  r  is  the  relative  degree  of  the  input/output  plant  dynamics. 
Finite-time-convergence  introduces  enhanced  robustness  to  perturbations  due  to  the  fact  that 
convergence  rate  in  the  vicinity  of  the  origin  will  be  higher  than  that  of  asymptotic  linear 
convergence. 

The  second  step  in  the  SMC  design  is  to  design  a  control  law  to  provide  convergence  to 
and  stability  on  the  sliding  surface  c  =  0 .  The  control  law  [8]  usually  consists  of  two  additive 
terms 

^  —  ^eq  ^un  •  (  ^  •  ^  ) 

In  case  of  full  information  of  the  plant  and  external  inputs,  the  term  ueq ,  called  equivalent 

control  [8],  can  be  determined  to  provide  exact  keeping  a  =  0  .  In  uncertain  case  it’s  substituted 
by  an  estimate,  ueq ,  of  any  portion  of  ueq ,  which  can  be  estimated  (calculated).  All  the  rest  of 
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uncertainty  is  captured  by  another  part,  uun ,  which  is  additive  and  complimentary  to  ueq  and 
can  capture  even  totally  uncertain  ueq ,  accepting  ueq  =  0 . 

For  conventional  sliding  mode  a  =  0  ,  the  examples  of  uun  are 

•  discontinuous  control  [8] 

um  =  />sgn(cr) ,  (1.2) 

•  continuous  control  [16] 

uun=p-lfcT.  (1.3) 

For  the  second-order  sliding  mode  &  =  a  -  0  continuous  control  law  [14]  is 

uun  =  Pi  +  Po  jsgn(o-)^ r.  (1.4) 

W  ' 

All  the  uncertainties  of  the  plant  are  accumulated  in  a  single  quantity,  the  relative  uncertainty  of 
equivalent  control  uecj  -  ueq  <  L ,  which  should  be  bounded  in  a  reasonable  domain  of  its 

arguments.  The  relative  uncertainty  gives  the  value  for  coefficients  p,  p,  P\ ,  Po  to  get  the 

sufficient  control  authority  to  provide  existence  of  sliding  mode.  The  greater  relative  uncertainty, 
the  more  control  authority  is  required  to  guarantee  compensation  for  it.  (The  method  to  get 
rigorous  proof  to  existence  of  sliding  mode  [8]  is  based  on  the  Lyapunov  function  technique.) 

1.3  SMC  Implementation  Issues 

In  practical  implementation  of  the  control  law  (1.2)  or  (1.4),  due  to  finite  time  switching 
of  sgn(-) ,  the  very  high  frequency  auto-oscillation  regime  is  established.  It  has  the  ability  to 
capture  all  the  uncertainties  and  disturbances  and  cancel  out  their  effect  on  the  closed-loop 
dynamics.  It  has  robustness  to  asynchronous  sample  rates  of  measurements  and  ZOH  control  in 
digital  controllers  as  well.  Control  (1.2)  provides  for  accuracy  of  a  -  0  to  be  proportional  to  r 
(switching  time  of  sgn(-)).  The  continuous  control  (1.4)  provides  for  accuracy  of  cr  =  0  to  be 

proportional  to  r~.  The  simplest  control  law  (discontinuous  (1.2)  or  continuous  (1.3))  provides 
for  given  robustness.  All  further  sophistication  (for  control  to  operate  at  almost  maximal 
amplitude  and  rate  without  saturation,  or  to  provide  the  output  feedback  only,  instead  of  state 
feedback)  increase  performance  and  retain  robustness. 


The  discontinuous  SMC  is  relevant  for  missile  autopilot 

There  are  two  major  objections  against  discontinuous  control  law:  the  on-off  switchable 
mode  is  not  energy  saving  and,  in  case  when  control  is  the  fin’s  deflection  angle,  they  cannot 
move  instantly.  However,  discontinuous  control  law  is  particularly  applicable  to  the  pulse  type 
on-off  thrust  rocket  engines  with  pulse  frequency,  which  is  much  higher  than  the  bandwidth  of 
the  overall  system  dynamics,  where  the  PWM  (Pulse  Width  Modulation)  or  PWPF  are  employed 
so  far  [7]  to  produce  input  equivalent  to  the  continuous  control  law.  Besides,  many  amplifiers  in 
electromechanical  actuators  work  in  PWM  mode  to  transmit  the  control  signal  to  the  amplified 
actual  power  signal.  Instead  of  PWM,  a  SMC  will  automatically  produce  high  frequency  auto¬ 
oscillation. 

The  control  structure  (1.1)  leaves  the  place  for  any  parallel  adaptive  corrections  (based  on 
explicit  or  implicit  on-line  identification).  The  less  ueq  becomes  uncertain,  the  better.  The  direct 

link  of  uncertainty  bounds  to  control  authority  gives  the  room  for  further  optimization  of 
specifications  on  control  effectors  design.  In  case  of  severe  failure  or  condition  change,  which  is 
above  the  range  specified  initially,  the  automatic  reconfiguration  can  be  made  [17]  without 
change  of  main  structure. 

At  last,  one  should  mention  that  SMC  design  is  less  theory  intensive  than  many  novel 
approaches,  such  that  a  customer  can  operate  independently  with  the  new  technology  being 
delivered.  Thus,  the  simple  fixed  controller  structure  for  a  plant  with  given  relative  degree  and 
the  given  range  for  cumulative  relative  uncertainty,  and  clear  distinct  steps  to  design  a  controller 
and  to  tune  it  up,  make  SMC  to  be  a  promising  technique  for  practical  missile  autopilot  design. 

In  this  report,  we  develop  a  multiple  loop  control  system  featuring  sliding  modes  in  each 
loop  and  demonstrate  the  application  of  this  SMC  design  to  a  generic  missile  autopilot 
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2.  Missile  Autopilot  General  Structure 

The  first  proposal  of  this  work  is  to  establish  a  new  autopilot  structure,  and  then  to  apply  all  the 
corrections,  adaptations,  gain  scheduling,  if  necessary.  For  decades,  a  classical  generalized  linear  autopilot 
structure  to  track  the  commands,  issued  by  a  guidance  law,  has  had  the  form  presented  in  Fig.  1  [3] 


Servo 


Fig.  1  Classical  Autopilot  Structure 

Where,  the  stability  loop  is  to  stabilize  the  servo  in  case  it’s  a  pure  integrator;  the  damping  loop 
is  necessary  to  restrict  oscillations  and  overshoots  to  provide  adequate  system  stability  margins; 
the  control  loop  is  responsible  for  command  signal  tracking. 

In  the  missiles,  the  electrical  control  signal  furnished  by  the  controller  is  superimposed  on 
a  direct  rate  gyro  feedback  implementing  damping  loop.  Sometimes,  a  linear  feedback  controller 
is  accomplished  by  adding  a  lead  compensator  to  the  forward  loop  to  stabilize  the  unstable 
airframe  [2],  Many  new  linear  and  nonlinear  design  approaches  have  left  the  first  two  loop 
without  change  and  redesign  only  the  third,  servo- tracking  control  loop  [1 1],  As  far  as  the  SMC 
enforces  the  closed-loop  dynamics  with  parameters  of  entirely  designer’s  choice  and  the  order 
being  equal  to  the  relative  degree  of  input/output  dynamics,  this  classical  structure  is  no  more 
necessary.  Moreover,  to  provide  the  extreme  agility  of  a  missile,  this  artificial  damping  loop  will 
be  neutralized  by  a  SMC  and  will  not  be  effective  anymore.  The  plant  can  even  have  an  unstable 
actuator,  since  all  the  problems  with  stability  and  performance  will  be  resolved  by  a  SMC. 

For  instance,  the  following  structure  for  attitude  angles  tracking  is  proposed  for  minimum 
phase  plants. 
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<p> 


C  ontr  1 


Contr  2 


-^0^]Contr  3 


Fig.2  SMC  based  multiple  loop  missile  autopilot  (General  Scheme) 

where  the  sliding  quantity,  cr ,  is  formed  in  each  loop,  and  enforced  via  virtual  control  that  is  to 
be  followed  by  in  the  next  inner  loop. 

In  the  next  section,  we  design  a  multiple  loop  SMC  system. 
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3.  Smooth  Sliding  Mode  Control  in  Multiple  Loop  Tracking  Systems 

3.1  Introduction 

In  multiple-loop  backstepping-type  control  systems,  an  important  issue  is  to  provide  a  so- 
called  virtual  control  signal  to  be  smooth  in  an  outer  loop,  for  it  has  to  be  tracked  by  inner 
cascades  of  a  multi-loop  system.  Smoothness  of  a  control  signal  is  not  easy  achievable  in  control 
systems  with  sliding  modes  without  degrading  high  precision  and  robustness. 

The  idea  of  this  work  is  to  use  combination  of  both  the  sliding  mode  estimator  for  a  plant 
disturbance  and  smooth  sliding  mode  controller  to  design  a  smooth  sliding  mode  controller  that 
is  robust  to  disturbances  and  provides  finite  time  convergence  to  the  custom  made  sliding 
surface. 

Two  types  of  smooth  sliding  mode  controllers  are  to  be  designed. 

•  First  order  smooth  finite  reaching  time  sliding  mode  controller  with  a  traditional 
sliding  mode  observer  for  disturbance  observation.  Under  the  discrete-time  control  with  a 
zero-order  hold,  the  accuracy  of  holding  the  trajectories  on  the  sliding  surface  is  of  the  first 
order  real  sliding  0(T ) . 

•  Second  order  smooth  finite  reaching  time  sliding  mode  controller  with  a  second 
order  sliding  mode  observer  for  disturbance  observation.  Under  the  discrete-time  control  with 
a  zero-order  hold,  the  accuracy  of  holding  the  trajectories  on  the  sliding  surface  is  of  the 
second  order  real  sliding  o(t2  ). 

For  many  control  applications.  Sliding  Mode  Control  (SMC)  [8-10]  has  been  proven 
efficient  technique  to  provide  high-fidelity  performance  in  different  control  problems  for 
nonlinear  systems  with  uncertainties  in  system  parameters  and  external  disturbances.  Ideal 
sliding  modes  feature  theoretically-infinite-frequency  switching,  while  the  real  conventional 
sliding  modes  feature  high  finite  frequency  switching  of  an  input  signal  (control).  Such  a  mode 
might  be  unacceptable  if  the  control  signal  has  to  be  tracked  by  inner  cascades  of  a  multi-loop 
system.  Trading  the  absolute  robustness  on  the  sliding  surface  for  the  system  convergence  to  a 
small  domain,  the  “boundary  layer”,  around  it  under  a  continuous  control  law,  the  methods  of 
this  group  employ  a  high-gain  saturation  function  [8],[  1 9]  or  [20]  a  sigmoid  function.  In  the 
continuous-time  control  systems  with  sampled-data  measurements  and/or  discrete-time  control 
action  (zero-order  hold  with  digital  control),  different  types  of  closed-loop  boundary-layer 
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dynamics  are  employed  to  provide  for  a  smooth  control,  varying  from  self-adaptive  saturation 
level  functions  [21]  to  fixed-gain  deadbeat  controls  with  disturbance  estimation  using  delayed- 
time  data  [22-24],  Another  alternative  [25]  to  the  latter  approach  is  to  incorporate  into  the 
“boundary-layer”  dynamics  an  exosystem  model  for  disturbances  [26]  avoiding  the  direct 
observer-based  disturbance  estimation. 

The  idea  to  hide  discontinuity  of  control  in  its  higher  derivatives  has  been  realized  using 
higher  order  sliding  modes  [14,27],  The  resulting  higher-order  sliding  mode  is  of  enhanced 
accuracy  and  robustness  to  disturbances.  However,  a  drawback  of  the  direct  application  of  this 
approach  to  chattering  attenuation  is  that  it  cannot  tolerate  unmodeled  fast  dynamics.  Therefore, 
the  designed  continuous  control  cannot  be,  for  instance,  an  outer-loop  feedback  in  a  multi-loop 
control  system. 

The  idea  of  this  paper  is  to  use  both  the  disturbance  estimation  and  the  higher-order 
sliding  mode  techniques  to  design  a  continuous  sliding  mode  control,  providing  finite-time 
convergence  to  the  sliding  surface  and  establish  the  second-order  sliding  mode  in  absence  of 
unknown  disturbances.  In  case  when  disturbances  are  present,  the  disturbance  observer 
determines  the  accuracy.  Employing  the  second-order  observer  [14],  the  second-order  sliding 
accuracy  can  be  achieved.  The  main  contribution  of  this  paper  is  in  further  development  of  the 
approach  presented  in  the  work  [28], 

3.2  Tracking  Problem  Formulation 

Consider  a  MIMO  plant  with  n  states  and  m  controls,  where  the  “diagonalization  method41 
[8]  has  been  applied  producing  m  independent  dynamics  for  each  input-output  channel.  Then, 
consider  the  following  SISO  nonlinear  uncertain  system  that  can  represent  any  input-output 
channel  (we  assume  relative  degree  is  equal  to  3,  although  the  given  approach  can  be  generalized 
to  rtb  order  system  as  well) 

i,  =  *2  +  tpx  (-Y,  )  +  /,(*,,/), 

x2  =x3  +  <p2(xl,x2)  +  f2(xl,x2,t),  ^  ^ 

x2=<p3(xx,x2,x3)  +  f3(-,t)  +  u, 
y  =  x  i 

where  the  <z>(  (•) -functions  are  known,  are  uncertain  time-varying  functions  that  are 

bounded  in  any  bounded  compact  set  of  their  arguments.  The  problem  is  to  provide  for  the 
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output  tracking:  y(t)  -»  yc(t ) .  The  design  problem  is  to  achieve  this  tracking  in  multiple  loop 

sliding  modes  using  full  state  feedback  and  backstepping  design. 

3.3  Multiple  Sliding  Surface  Design  with  lst-Order  Sliding  Modes 
A  multiple  sliding  surface  implementation  of  the  back-stepping  approach  looks  as  follows  [30], 

1)  Define  the  first  sliding  surface 

51  =  yc  -y  =  xf  -.r,  =  0,  (3.2) 

and  the  following  desired  closed-loop  dynamics  for  the  sliding  quantity  5, 

S|  =  _/i(Si)  +  S2  +  2j  (3.3) 

where  S2  =  x2 -x2 ,  x2(t)  is  to  be  defined,  Ij  =  (if  -/i(-,o)_(*f  -/i)>  xi  -  Tc»  and  /[(•)  is  such  a 
function  that  the  homogeneous  part  of  the  system  (3.3)  is  a  finite  time  convergent  equation. 

From  (3.1)-(3.3)  one  can  formally  obtain 

Sj  —  yc  ~ x2  —  <P\ (•)  —  f[ (•,/)  =  —  Y\ (Sj )  +  x2  —  x2  +  (ij  —  /, (-,o)- (-h  ~  f\)  (3.4) 

From  (3 .4)  we  derive  the  following  reference  command  for  the  next  step 

x2  =  if  -f\~<P\  (v,  )  +  /i(S,)  (3.5) 

where  if  ,/j  are  the  best  estimates  of  the  reference  signal  and  uncertainty  in  the  first  loop  that  are 

based  on  dynamic  observers. 

2)  The  second  sliding  surface  is  defined  as 

52  =  x2  —x2  =  0 ,  (3-6) 

Assuming  existence  of  (3.6)  and  finite  time  convergence  of  (3.3)  even  with  Ij  *  0,  the  sliding 
mode  on  the  surface  (3.2)  will  be  achieved. 

Sliding  quantity  S2  is  calculated  using  feedback  on  x2  and  equation  (3.5)  for  x2.  An 
alternative  way  to  produce  S2  using  S|  feedback  only  can  be  obtained  as  follows.  Derive  the 
following  identity  from  (3.1),  (3.2) 

X\  -  f\  -  S|  +  x2  +  (.V, ) , 

then  one  can  estimate 

X\  ~  f\  =  5]  +  x2  +  q>\  (.Yj ) .  (3.7) 

From  (3. 5), (3. 7)  we  obtain 

S2=S, +/,(£,).  (3.8) 
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There  are  many  ways  to  obtain  (3.8).  One  can  estimate  the  derivative  and  calculate  the  nonlinear 
term,  one  can  try  to  estimate  them  both  at  once;  at  last,  one  can  apply  the  nonlinear  DSM 
approach  and  enforce  (3.8)  in  the  system  motion  in  an  auxiliary  DSM. 

The  closed-loop  dynamics  for  the  sliding  quantity  S2  is  selected  in  the  form 


52  =  -r2(S2)  +  S3+?2  (3.9) 

where  S3  =  x$-x3,  xc3(t)  is  to  be  defined,  e2  =  (xf  -  f2{-,t))-[xc2  -/2),  and  y2 0  is  of  the  same 
class  as  (•) .  Similar  to  Eq.(3.5)  the  reference  command  xc2{t)  is  obtained 

*3  =  ^2  -  A  -  ^2 (-^1  ,-^2 )  +  ^2(^2)  -  (3-10) 

3)  Finally,  the  third  sliding  surface, 

53  =*3 -*3=0,  (3.11) 

is  achieved  under  the  control 

u  =  -<p3  (•)  +  *!, 3S3  +  /f0>3sgn(53) .  (3.12) 

The  closed-loop  system  in  the  ( SUS2,S3 )  state  space  is  derived  as 

51  =  ~/i(Si)  +  S2  +2j, 

52  =  ~/2(S2)  +  S3  +e2,  (3.13) 

53  =  x2  —  —  k  1 3S3  — 3  sgn(53). 


Dynamics  (3.13)  plus  the  error  dynamics  in  the  dynamic  observers  that  govern  (e,,e2)  complete 
the  total  closed-loop  dynamics  of  the  plant. 

In  the  work  [30],  it’s  proved  the  closed-loop  stability  and  convergence  of  (3.13)  with  first-order 
linear  dynamic  observers.  These  observers  have  served  as  integral  filters  to  smooth  out  possibly 

non-Lipshitz  behavior  of  if,  if,  *3  •  That’s  why  it  was  possible  to  define  feedback  terms  in  each 

compensated  dynamics  as  ^.(S,.)  =  kiAS,  +ki0sgn(Si),i  =  1,3 .  These  functions  provide  for  finite 

time  convergence  in  (3.13),  but  they  are  not  differentiable.  Linear  smoothing  filters  of  the  work 
[30]  overcome  the  problem  of  “explosion  of  terms”  trading  robustness  for  stability.  However,  in 
numerical  implementations  this  discontinuous  form  is  rarely  used,  for  even  with  filtering  the 
term  sgn(S)  behavior  in  outer  loops  is  not  good  for  inner  tracking  loops.  There  is  no  finite  time 
convergence  and  the  true  sliding  mode  [20]  in  this  case.  Another  approach  is  developed  in  this 
work  to  establish  finite  time  convergence  in  each  loop  to  a  dynamic  integral-type  sliding  surface. 
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Additionally,  the  integral  part  specifies  time  scale  separation  [17]  between  the  loops  to  ensure 
overall  stability  of  a  backstepping-type  tracking  system. 

3.4  Multiple  Dynamic  Sliding  Surface  Design  with  2nd-Order  Sliding  Modes 


3.4.1  Continuous  Control  for  Finite  Reaching  Time  Sliding  Mode 

Finite-reaching-time  continuous  standard-sliding-mode  controllers  have  been  developed 
in  works  [16], [28]  for  the  first  order  a  -dynamics 

&  =  f(a,t)  +  u .  (3.14) 

They  provide  for  finite-time-convergence  of  the  first-order  closed-loop  <x  -dynamics.  One  of  the 


forms  in  the  work  [  1 6]  is  given  by 

<7 


CT+  p- 


H 


0.5 


=  0, 


(3.15) 


see  Fig.3. 


<7 


Fig.3  Nonlinear  Terminal  Sliding  Manifold 
In  absence  of  uncertainty  in  the  function  / (a ,t) ,  the  control  law 

u(a)  =  -f(cr,t)~  (3.16) 

M  ' 

renders  the  closed-loop  dynamics  (3.15),  as  a  finite  time  convergent  nonlinear  manifold.  When 
the  function  f(cr,t)  is  totally  uncertain,  the  continuous  control  law 

"<<T)=-^r5r- 

M  (3.17) 

provides  for  convergence  to  the  arbitrarily  small  domain  of  attraction,  the  boundary  layer,  around 
the  sliding  surface  a  =  0  in  a  standard  sliding  mode,  where  the  gain  p  and  the  uncertainty  limit 

L  determine  the  boundary  layer  thickness.  The  drawbacks  of  this  controller  are  that  the 
uncertainty  limit  defines  the  boundary  layer,  and  even  in  absence  of  uncertainty  the  domain  of 
attraction  to  a  =  0  is  proportional  to  the  discrete  interval  T  under  the  discrete-time  control  (first- 
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order  sliding  accuracy).  An  additional  problem  for  a  mupltiple  loop  system  is  that  (3.17)  is  not 
smooth  enough  to  be  r  times  differentiable. 

3.4.2  Conditions  on  Smoothness  of  a  Virtual  Control 

In  order  for  the  control  of  form  (3.17)  to  be  a  virtual  control  to  be  followed  by  inner 
cascades  of  total  order  r  of  a  multi-loop  system,  it  has  to  be  r  times  continuously  differentiable. 

In  work  [16],  finite  time  convergence  has  been  proved  for  the  closed-loop  dynamics 

<r  +  /?|cr|a  sgn(cr)  =  0 ,  (3.18) 

where  a  e  (0,l).  Similar  to  (3.17),  it  gives  us  the  virtual  control  in  the  form 

U(\)=-p\o\a  5gn{a) .  (3.19) 

For  to  be  followed  by  a  first  order  tracking  system  is  has  to  be  one  time  continuously 

differentiable.  From  (3. 18), (3. 19)  we  have 

«(!)  =  -p2a\o\2a  1  sgn(cr) .  (3.20) 

Eq.  (3.20)  gives  the  following  condition  on  smoothness  of  the  virtual  control  at  the  origin 
a>l- 

In  case  of  an  rlh  order  tracking  system,  we  have 

u(r\r )  ~  -\o\{r+l)a~r  sgn(er) .  (3.21) 

Finally,  combining  conditions  of  finite  time  convergence  of  the  compensated  dynamics  and 
smoothness  of  the  virtual  control,  we  obtain  the  following  condition  for  the  rth  order  tracking 
system. 

Condition  on  Smoothness.  Terminal  sliding  mode  can  be  enforced  on  the  surface  (3.18)  by  an 
r,h  order  tracking  system  in  all  loops  if 

—  <a<\.  (3.22) 

r  + 1 

Applying  this  condition  to  system  (3.1)  we  must  have  ycarcl  if  we  select 
/i(5i)  =  |.S’1|<*  sgnC^i)  for  (3.3).  One  should  not  worry  about  overall  stability  of  system  (3.13) 
and  convergence  rate  of  the  estimation  error  (ij,^) -dynamics  if  time  scale  is  applied  in  each 
loop.  Multiple  time  scale  is  achieved  in  dynamic  sliding  surfaces  designed  in  the  next  section. 
Second  order  sliding  modes  are  established  on  dynamic  sliding  surfaces  that  govern  compensated 
time  scaled  dynamics  in  each  loop. 
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3.4.3  Multiple  Time  Scale  Dynamic  Sliding  Surfaces 
The  first  (outer)  loop  design 

Consider  the  tracking  problem  for  system  (3.1).  We  introduce  the  following  dynamic 

sliding  surface  for  the  compensated  dynamics  in  the  outer  loop  of  the  expected  3-loop  system 

<J\  +  px  J|cr, sgn(cr,  )dt  =  el+cl  j\e{\ai  sgn(e, )dt ,  (3.23) 

where  ex  is  the  input  and  the  dynamic  sliding  surface  quantity  ax  is  the  output.  If  the  sliding 
mode  exists  on  the  surface  <jx  -  0 ,  then  from  (3.23)  we  obtain  in  sliding  mode 
ex  +  Cj  j\ex\ai  sgn(e{ )dt  =  const ,  ax  =  0 , 


or 

ex  =-Cx\ex\a]  sgn(ej),  (3.24) 

which  is  a  finite  time  convergent  system,  where  the  terminal  time  is  a  function  of  ex(0 ),ax,cx . 

To  enforce  convergence  to  =  0 ,  we  consider  <7,  -dynamics  using  Eq.(3.1)  and  introducing  the 
virtual  control  x2 ,  and  the  tracking  error  e2  =  x2  -  x2  of  the  inner  loop 

<7 1  +/?1|cr1|^J  sgn((7j)  =  yc  - <P\{x{) -  - (x2 -e2)  +  cx \ex\ai  sgn(e, ) .  (3.25) 

The  virtual  control  x2  is  designed  as 

x2  =  h-<P\(x\)- f\i-d)  +  cx\ex  T1  sgn(ej)  +  />0  J[o"  j  |^’°  sgn(<7j)Jt ,  (3.26) 

Given  (3.26)  and  estimation  error  ex  =  (yc  -  fx  (.,/))  -  (yc  -  f\ (.,/)) ,  the  outer  loop  cr,  -dynamics 
closed  under  control  (3.26)  is  obtained 

+  P\  | <7 x  l^'1  sgn(<7, )  +  Pq  J|<7,  |^  °  sgn(<7[  )dt  =  ex+e2.  (3.27) 

Assuming  exact  estimation,  i.e.,  ex  =  0,  and  fast  convergence  of  e2  to  zero  (time  scale  of  the 
inner  loop  dynamics),  we  have  the  second  order  cq  -dynamics 

<X\  +  P\ ,  -rrr1 - +  A  ki \A’°  sgn(<7, )  =  0 .  (3-28) 

sgn(^) 

which  is  a  finite  time  convergent  differential  equation  given  conditions  formulated  in  the 
following  lemma. 

Lemma  1  Consider  the  following  <j -dynamics 


&  =  -ax\of ' 2  sgn(cr)  -a0j\o j1 7  3  sgn(tr)^  r , 


(3.29) 


ax  >  0,^0  >  0 ,  which  can  be  equivalently  presented  by  the  system  of  two  first-order  equations 
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(3.30) 


*1  ~  x2~a\\x§ 12 sgn(Ar, ), 

*2 =-^okl1/3sgn(^i). 

where  xx  =  a.  The  system  (3.29)  is  asymptotically  stable  and,  moreover,  it  approaches  the  origin 
in  a  finite  time. 

Proof:  It’s  not  difficult  to  prove  asymptotic  stability  of  system  (3.29),  let  a  Liapunov  function 
candidate  be 

2  -vl 

V(xx,x2)  =  Y+  lao\z\ U2sgn(z)dz, 
o 

V(x)  >0,  if  *  e  \  {0} ,  then  the  Liapunov  function  derivative  will  be 

~  1/2 

•  dv  x2-ax\xx\  s§n(*l)  l  |L4  A  .p  ~2  wm 

V- -  =  -&oal\x\r  2<0,iixe9?  \ {0}  . 

dx  _-a0|xi|,/3sgn(jfi)  J 

However,  in  order  to  prove  the  finite  time  convergence  to  the  origin  of  system  (3.29),  one  have 
to  apply  special  topics  of  Liapunov  analysis  of  the  finite  time  convergence  differential  equations. 
In  the  work  [29]  the  following  theorem  has  been  proved. 

Theorem  Let  xe  Da  0?" ,  .i-  =  f(x) ,  / : -*  dV  is  continuous  on  an  open  neighborhood D 
of  the  origin  and  locally  Lipshitz  on  Z>  \  {0}  and  /(0)  =  0.  Suppose  there  exists  a  continuous 
function  V :  D  -» such  that  the  following  conditions  hold 

(i)  V  is  positive  definite; 

(ii)  V  is  negative  on  Z>  \  {0} ; 

(iii)  there  exist  real  numbers  k  >  0  and  a  e  (0,1) ,  and  a  neighborhood  N  a  D  of  the 
origin  such  that  V  +  kVa  <0  on  N\{0} . 

Then,  the  origin  is  a  finite-time-stable  equilibrium  of  x  =  f(x) . 
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The  function  V  is  satisfying  (i),(ii)  of  the  theorem.  To  prove  (iii)  we  consider  the  validity  of 


(  J2 


V  +  kVa  =  -<2'0tf'i|*i|6  +k 


y+i^oFi 


a 


<0. 


The  inequality  can  be  transformed  to 


'2  ^  \  1  5  | 

/ 


(3.31) 


In  a  small  neighborhood  of  the  origin,  where  x2 « l,*i « 1,  and  xi-»0,  i.e. 
x2  —  l^i |* 7 2 sgn(jcj )  ->  0,  we  have  |x2|  ~  |^i|I/2 ,  then  from  (3.31)  we  conclude  that  the  left  part  in 
(3.3 1)  will  be  dominated  by  |xj| ,  when  x2  « l,*i  « 1 ,  while  the  right  part  will  be  proportional  to 

|x]|6ff.  There  exists  <zg(0,1)  such  that  — —  <1,  hence  one  can  always  select  k>  0  and  a 

6a 

1  5 1 

neighborhood  TV  c  D  of  the  origin  such  that  (aQai)a\xi\^a  will  be  dominated  over  the  right  part 
in  (1 1)  which  is  proportional  to  |xj|  in  a  small  N  c  D .  Thus,  the  proposition  (iii)  of  the  theorem 
is  true,  and  system  (3.29)  has  a  finite-time  stable  equilibrium  at  the  origin  for  any  ax  >  0,ar0  >  ^ 

Thus,  if  ,  =l/2,/?10  =1/3,  system  (3.28)  converges  to  the  origin  in  a  finite  time 
featuring  the  second  order  sliding  mode  di  =  &i  -  0.  In  the  sliding  mode,  the  compensated 

tracking  error  dynamics  is  governed  by  Eq.(3.24)  under  perfect  tracking  x2  — >■  x2  . 

To  avoid  explosion  of  terms  the  virtual  control  (3.26),  i.e. 

x2  =  >'c  -  <P\  (*1 )  “  fi 0,0  +  q  \el  I"1  sgn(e, )  +  p0  J| crx \A-°  sgn(cr,  )dt , 
must  be  continuously  differentiable  twice.  This  gives  us  design  condition  on  at 

§  <  <  1 , 

and  the  following  condition  must  be  checked  for  .  For  the  second  derivative  of  x2  be 
bounded,  i.e.  since  x%  ~ - — ,  we  must  have 


^1 


pi 


<  00  .  From  (3.28)  we  obtain  that 
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The  second  loop  design 

Dynamics  of  the  tracking  error  e2  =  x2  -  x2  is  to  be  enforced  in  the  second  loop  on  the 
dynamic  sliding  surface 

a2+ pi  j\o-2f2']  sgn(cr2)dl  =  e2+c2  j\e2 \a2  sgn (e2)dt ,  (3.32) 

which  is  similar  to  the  first  one  (3.23)  with  one  difference,  it  must  be  faster  enough  to  enforce 
sufficient  time  scale  to  consider  e2  =  0  in  the  outer  loop. 

To  enforce  convergence  to  a2  =  0,  we  consider  cr2 -dynamics  using  Eq.(3.1)  and  introducing 
the  virtual  control  x2 ,  and  the  tracking  error  e2  =x2  -  x2  of  the  inner  loop 

&2  +/?2,i|cr2|^'1  sgn(<x2)  =  jf  -<p2{x\,x2)~  f2{.,t)-{x2  ~e3)  +  c2 \e2\ai  sgn(e2). 

(3.33) 

The  virtual  control  x2  is  designed  as 

*3  =  xc2-<p2{x{,x2)- f2{.,t)  +  c2\e2\a2  sgn (e2)  +  p2  0  J|cr2|^-0  sgn{a2)dt ,  (3.34) 

Given  (3.33)  and  estimation  error  e2  =(x2  -  f2 (.,/)) ~{x2  -  f2 (.,/)),  the  second  loop  a2- 
dynamics  closed  under  control  (3.34)  is  obtained 

&2  +  P2,\  sgn(^2 )  +  p20  j\cr2f2’0  sgn(cr2  )dt  =  e2+e2.  (3.35) 

Assuming  exact  estimation,  i.e.,  e2  =  0 ,  and  fast  convergence  of  e3  to  zero  (time  scale  of  the 
inner  loop  dynamics),  we  have  the  second  order  finite  time  convergent  cr2  -dynamics. 

The  third  loop  design 

The  very  inner  loop  is  to  be  designed  for  actual  control  u  to  stabilize  the  tracking  error 
e3  -dynamics 

"  <Pi ) -  fi (.,/) - u  .  (3.36) 

The  sliding  mode  control  is  designed  in  a  standard  way  [8] 
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(3.37) 


u  =  <p3(x],x2,x3)  +  psgn(e3) ,  p> 


*3 


+ 


Additional  requirements  must  exist  for  a2  in  order  for  if  to  be  bounded.  From  the  obtained 
condition  on  smoothness  we  have  ^  <  a2  <  1 . 

The  closed  loop  dynamics  in  the  third  loop, 

<?3  =*3C  -/3(->0-/?sgn(e3).  (3.38) 

is  a  finite  time  convergent  system. 


3.4.4  Conclusions  on  Multiple  Time  Scale  Dynamic  SMC  Design 
A  three-loop  tracking  control  system  has  been  designed  for  a  third  order  uncertain  SISO 
system.  Tracking  error  dynamics  is  enforced  in  each  loop  in  terminal  dynamic  sliding  surfaces. 
The  second  order  sliding  performance  is  provided  for  sliding  modes  on  dynamic  sliding  surfaces 
in  the  outer  and  the  inner  loops  under  virtual  controls.  Time  scale  separation  ensures  overall 
stability  of  Eqs.(3.27).(3.35),(3.38).  When  ex  =e2  =  0,  exact  tracking  y  -  yc  is  achieved  in  a 

finite  time  by  system  (3.1)  under  control  (3.37). 

A  signal  differentiator  and  a  disturbance  observer  must  accompany  the  presented  design 
to  ensure  e{  0,<?2  0  as  close  as  possible.  The  design  of  a  signal  differentiator  and 

disturbance  observer  on  its  basis  is  presented  in  the  next  section  on  the  example  of  calculating 

h  and  f\{.,t)  ■ 
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4.  An  Auxiliary  Tool:  Robust  SMC  Differentiator 
2-Sliding  Exact  Differentiator 


The  2-sliding  exact  differentitator  developed  in  the  work  [14]  can  be  employed  to 
produce  the  time  derivatives  of  a  smooth  reference  signal.  An  exact  differentiator  of  a  given 
signal,  x(t) ,  has  the  form 


Xy  =  Jx2dr, 

<  «  (X(t)-Xy)  r  ,  ,  ,  .  v  ,  (4.1) 

x2  =  P\ : - yy  +  Poj sgn{x(t)  -xx]dz. 

|x(0-x,|  ' 

where  the  output  x2  converges  to  x(t)  in  a  finite  time.  To  verify  this  statement,  we  can  identify 
the  dynamics  of  estimation  error,  e  =  x  -  x ,  using  Eq.  ( 1 6)  as 


e  +  ^  +  a  sgn(<?)  =  x(t).  (4.2) 

2  M 

This  is  a  finite-time-convergence  second-order  differential  equation,  which  provides  for 
2-sliding  mode,  e  =  e  =  0,  in  case  when  the  parameters  a,  A  are  appropriately  selected  in 

accordance  with  the  upper  boundary  for  |x(t)|  <  L .  A  proof  to  this  statement  gives  the  following 

lemma,  which  summarizes  the  results  of  the  work  [14]  about  2-sliding  modes  with  finite  time 
convergence. 


Lemma  2  Given  2nd  order  nonlinear  differential  equation  (4.2)  and  conditions  |.V(/)|  <  L  ,  a>4L, 

A  >  0.5Z.0'5  ,  e(0)  =  0,  e(0)=zo,  z'0  =  arbitrary’  constant,  any  its  solution  converges  to  the  origin  in 
finite  time. 

Remark  The  observer  (4.1)  is  subjected  to  the  chattering  of  the  output  for  noisy  signals. 
Therefore,  in  practical  observer  of  this  form,  a  low-pass  filtering  of  x2  is  required  for  a  smooth 
estimate  of  x(t) ,  as  shown  in  the  work  [14]. 

2-SIiding  Disturbance  Observer 

We  consider  disturbance  estimation  on  the  example  of  the  first  equation  in  (3.1),  i.e., 
fy  (.,/) .  The  state  observer  for  the  system 

Xy  =  <Py(Xy)  +  X2  +/,(.),  (4.3) 


23 


is  designed  as 


*1  =  <P\  (*1 )  +  *2  +  Pi  M° 5  s§n(e)  +  Po  jsgn (e)dt,  (4.4) 

where  e-xx  -xx .  The  observer  error  dynamics  is  derived 

e  =  <pl(xx)~  (px (x, )  +  /, (.)  —  u0,  uQ  =  P\ |e|° 5  sgn(e)  +  /?0  \sgn{e)dt (4.5) 
This  is  a  finite  time  convergent  dynamics  given  the  restriction  \<px  (^!  )-<px  (x!  )||  <  C||e||  for  all 
times.  Thus,  uQ  is  an  estimate  of  /j .  To  provide  smooth  derivative,  the  estimate  u0  is  low-pass 
filtered  before  entering  the  virtual  control  as  fx 


5.  Simulation  Example 


To  illustrate  the  disturbance  cancellation  characteristics  of  the  developed  method,  we 
consider  the  simplified  model  of  a  ballistic  interceptor  missile.  A  simplified  model  in  the  pitch 
plane  is  given  by 


v\  =-—(«„  +  u<l)sm9  +  — 

m  m  ,  (5.1) 

v2  =—  (ua+ud)cos9  +  ^--g 
m  m 


where  V  =  JvJ  +  v.2  is  the  velocity  vector,  m  is  the  mass  of  the  interceptor  in  kg,  y/x  and  y/2  are 
disturbances,  g=9.81m/s2  is  gravity,  9  is  the  pitch  angle  of  the  interceptor,  ua  is  the  attitude 
control  and  uj  is  the  divert  control. 

A  simplified  rotational  motion  model  for  the  interceptor  is  given  by 


9  =  q, 


(5.2) 


2 

where  q  is  the  pitch  rate  (rad/s),  J  is  the  moment  of  inertia  of  the  interceptor  (kgm  ),  r  is  the 
distance  between  application  of  the  attitude  control  force  ua  at  the  interceptor’s  center  of  gravity, 
and  y/m  is  the  rotational  disturbance  term.  The  output  equations  are  given  by 


(y,  =  4,  =rv 

}y2  =a  =  e-y' 


(5.3) 


where  y  is  the  missile’s  flight  path  angle,  a  is  the  angle  of  attack,  and  A„  is  the  normal 
acceleration.  Our  goal  is  to  design  the  attitude  control  thrust  to  stabilize  angle  of  attack  to  zero, 
assuming  that  the  divert  control  thrust  delivers  desired  normal  acceleration. 

lim|y,//)-y/tj|  =  lim|4r -A„\  =  0 

I  — >oo 1  /— *oo  / r  £ \ 

lim  I  \\r(t)  -  y2  (t)\  =  lim|0  -  a\  =  0 

where  Anr  is  a  desired  reference  profile  for  normalized  acceleration.  The  reference  profile  for  the 
normalized  acceleration  is  generated  by  the  guidance  system  to  minimize  miss  distance  in  the 
interceptor  using  a  number  of  methods,  such  as  proportional  navigation,  augmented  proportional 
navigation  and  other  advanced  guidance  laws  [1].  Forcing  o' to  zero  aligns  the  missile  body  with 
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the  flight  path  angle.  We  consider  the  normal  aerodynamic  lift  force  and  pitching  moment  as 
disturbances  entering  the  plant 

p  v2  p  V2 

N  =  ^—S-Cn  -a,  M  =  — — Sl-Cm  a, 

2  2  a 

y/x  =  -TV  sin/,  y/.  -  Ncosy,  y/a  -  M  . 

Finally  our  numerical  model  is  obtained  as 

Vx  =—((“«  +1000)  sin 

vz  =  ^((Ma  +  1000)cos<9+  {^z)-9.81, 

0  =  q, 

q  =  -0.\(ua+y/a), 


a  =  0-y. 


y  =  tan 


-l 


A 


\vx  J 


(5.5) 


v  =  i 


2  2 
v;  +  v;. 


y/x  =  -3-\0~4V2asiny,  y/z  =  3-l0~4V2acosy,  y/a  = -0.5 -\0~°V^a, 
vx (0)  =  2000m / s,  v. (0)  =  0 m/s,  0(0)  =  0.3 rad,  q( 0)  =  0 rad /s. 

The  plant  output  to  be  stabilized  is  angle  of  attack,  the  attitude  control  trust  is  actuated 
via  a  first  order  actuator 

ua=-20(ua-u),  (5.6) 

so  the  actual  control  signal  is  u,  and  we  have  a  third  order  input-output  dynamics 

A„ 


6  t/2  . 


a  =  q - , 

V 


(5.7) 


q  =  -0.\(ua+y/a), 
ua  =  -20(wa  -  u), 

where  we  consider  An  =  -vx  sin/  +  v.  cosy  and  y/a  as  disturbances  in  the  first  and  second 

loops  of  a  3-loop  control  system  to  be  designed. 

The  stabilizing  controller  is  designed  as  follows.  The  first  dynamic  sliding  surface  is 
selected  as 
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(5.8) 


Cj  +0.5  J| ctj | 1 7 2  sgn((Ti)dt  =  +- 0.05  J|<3r|2/3  sgn(a)dt  . 

Then,  <rx  -dynamics  are  identified,  introducing  eq  =  qc  -  q , 

61  =  — O-SjcT!  |1/2  sgn(oT! )  +  0.05|<ar|2/3  sgn(a)  +  qc  -eq  --y. 

Virtual  control  in  the  first  loop  is  designed 

qc  =  -0.05|ff|2/3sgn(«)-1.5j|<7,|1/3sgn(<x1)<fr .  (5.9) 

The  second  dynamic  sliding  surface  with  appropriate  time  scale  is  selected  as 

1 1  /  2 

cr2 +3j|cr2|1/2sgn(cr2)^  =  e9 +0.5  jje?|  sgn  (eq)dt.  (5.10) 

Then,  cr2  -dynamics  are  identified,  introducing  eu  =  uac  -  ua , 

1/2 

62  =  -3|ct2|1/2  sgn(cr2)  +  0.5e9  sgn{eq)  +  qc  +0A(uac -eu)-0Ai//a. 

Virtual  control  in  the  second  loop  is  designed 

uac  =  10  -\-qc  -0.5|e?|1/2  sgn(eg)  -10  j]cr2|1/3  sgn(<r2)<*j  .  (5.1 1) 

Tracking  error  dynamics  in  the  third  loop  are  identified 
eu  -  uac  +  20 ua  -  20 u  . 

Actual  discontinuous  control  is  designed 

u  =  20sgn(e„ ) .  (5.12) 

Finite  time  convergent  multiple  sliding  surfaces  dynamics  are  obtained 

61  +  0.5|cr,|'/2  sgn(cr,)  +  1.5  J|cr,[1/3  sgn{<7x)dt  =  -eq  ~y~, 

62  +3|cr2|l/2  sgn(cr2)  +  10j|cr2|1/3  sgn (a2)dt  =  -eu  -0.\y/a,  (5.13) 

eu  =  uac  +  20 ua  -  400sgn(ei( ). 

Analyzing  (5.13),  one  can  conclude  that  eu  reaches  zero  in  a  finite  time  robustly  to  bounded 
behavior  of  iiac ,  ua  .  When  eu  =0,  cr2  approaches  in  a  finite  time  a  small  domain  around  zero 
attenuating  the  disturbance  i//a .  When  <x2  a  0 ,  eq  goes  to  zero  in  a  finite  time  according  to 

eq  =-0.5|e?  °'5sgn(e(?).  When  eq  ~0,  crl  approaches  in  a  finite  time  a  small  domain  around 
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zero  attenuating  the  disturbance  ~ .  When  erx  »  0 ,  a  goes  to  zero  in  a  finite  time  according  to 
a  =  -0.05|cr|2/3  sgn(«r) .  Results  of  a  simulation  are  given  in  Figs. 
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6.  Scientific  Merit  and  Impact  of  the  Research 

•  The  new  smooth  sliding  mode  controllers  ihat  employ  first  order  and  second  order 
sliding  modes  and  sliding  mode  observers  are  proposed  to  stabilize  multiple  sliding 
surfaces  in  the  sliding  mode  controlled  plant. 

•  The  proposed  smooth  sliding  mode  controllers  can  be  used  is  control  systems  where 
control  actions  must  be  smooth,  in  particular,  in  multiple  loop  control  systems. 

•  Under  the  assumption  that  the  observer  captures  the  uncertainty  and  cancels  its  effect 
on  the  closed-loop  dynamics,  the  system  motion  in  the  vicinity  of  the  sliding  surface 
is  governed  by  a  second-order  finite  time  convergent  dynamics  that  establishes  a 
second-order  sliding  on  the  sliding  surface.  Thus,  both  parts  of  the  control  law 
(feedback  Pi-type  control  and  disturbance  observer)  are  the  second  order  SMC’s. 
Under  the  discrete-time  control  with  a  zero-order  hold,  the  accuracy  of  holding  the 

trajectories  on  the  sliding  surface  is  of  the  second-order  real  sliding,  o{t2  ). 
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